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Abstract: We point out that quintessence with an exponential potential 
y oe -/3Q/v3M p can accoun t f or the present observed acceleration of the universe, with- 
out necessarily leading to eternal acceleration. This occurs for 2.4 < (3 < 2.8. Thus 
a cosmological horizon, which is supposed to be problematic within the context of 
string theory, can be avoided. We argue that this class of models is not particularly 
fine-tuned. We further examine this question in the context of a modified Fried- 
mann equation, H 2 oc p + p, which is suggested by higher dimensional self-tuning 
approaches to the cosmological constant problem. It is shown that the self-tuning 
case can also be consistent with observations, if 1.8 < (3 < 2.4. Future observations 
of high- z supernovae will be able to test whether (3 lies in the desired range. 
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1. Introduction 

There is convincing evidence that the universe is presently dominated by a form of 
dark energy density which is decreasing significantly more slowly with time than the 
energy density of ordinary matter. The first Doppler peak in the cosmic microwave 
background fluctuations is strongly consistent with a flat universe whose density 
is the critical one p c , yet the matter energy density is known by other means to 
be no more than Q m = p m J p c ~ 0.3 0. Recent observations of the higher peaks 
by the various CMB experiments give Q m = 0.25 ||. The Hubble diagram deduced 
from high redshift type I supernovae provides independent evidence for an additional 
component f^A of the energy density Q. From an empirical viewpoint, pure vacuum 
energy (A) is the simplest explanation, but theoretically it is difficult to explain why 
A is some 124 orders of magnitude smaller than the natural scale set by the Planck 
mass, Mp. If the dark energy is due to a rolling scalar field Q, quintessence [|5], whose 
potential energy vanishes as Q — > oo, this might be a more natural explanation for 
why the dark energy density is small. 

A further motivation for quintessence could be coming from string theory, be- 
cause of its apparent incompatibility with de Sitter space || . An eternally accelerat- 
ing universe, which would result from a positive cosmological constant, seems to be 
at odds with string theory, because of the impossibility of formulating the S-matrix. 
In de Sitter space the presence of an event horizon, signifying causally disconnected 
regions of space, implies the absence of asymptotic particle states which are needed 
to define transition amplitudes. Quintessence, on the other hand, would seem to 
offer the possibility of temporary acceleration to account for current observations, 
without necessarily making the scale factor of the universe accelerate forever. 
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Recent papers have pointed out that, in fact, quintessence generically does lead to 
eternal acceleration, if the universe is accelerating now [0, |j (see also 0). Therefore 
it is concluded that string theory is equally at odds with quintessence or a positive 
cosmo logical constant. The purpose of the present paper is to explore some simple 
loopholes to this conclusion. 1 

Our initial motivation was an interesting possibility coming from the brane- 
world scenario, in particular, attempts to address the cosmological constant problem 
through self-tuning solutions to the Einstein equations flIT |. In this approach, it is 



assumed that our universe is a 3-brane with arbitrary tension A, embedded in an 
extra dimension. A scalar field living in the extra dimension adjusts itself so as to 
yield a static solution to the Einstein equations, regardless of the value of A, which 
otherwise would act like the 4-D cosmological constant and lead to inflation of the 
brane. Although there are many problems with this idea ]12), ref. [13| explored 



the question of how cosmology would be affected for brane observers assuming an 
acceptable model of self-tuning was found. If the scalar couples only to the volume 
element y r gd A x of the three- brane, [|13[ showed that it is possible to obtain a modified 
Friedmann equation of the form 

H 2 =(-) 2 = 2nG(p + p) (1.1) 



plus corrections of order G(p + p) 2 /TeV 4 (assuming the 5-D quantum gravity scale 
is of order TeV). This is a very interesting twist on normal cosmology because (1) 
it leads to no expansion in the case of vacuum energy, where p = —p\ (2) it is 
indistinguishable from the normal equation during the radiation dominated era, since 



p = p/3; and (3) it gives a Hubble rate only y 3/4 = 0.87 times smaller than normal 
during the matter dominated era, which would be difficult to distinguish from the 
standard value given the uncertainties on Q m (the fraction of the critical density in 
matter) and Ho (the present value of the Hubble parameter). Since the modified 
Friedmann equation eliminates conventional inflation, it would remove the obstacle 
to defining an S-matrix in the 4-D universe when A > 0. We might also expect it to 
ameliorate the problem with eternal acceleration in quintessence models. 

Moreover, if the equation of state for the dark energy, w — p/p, turns out to be 
w > — 1, as is still allowed by the supernova data, this could be indirect evidence for 
a modified Friedmann equation. Taking the time derivative of (|1 . 1| ) , one can show 
that the acceleration is 

~ = -TrGp(l + w)(l + Zw) (1.2) 
a 

so that the condition for positive acceleration is —1 < w < —1/3. This is to be 
contrasted to the standard result, - = — (47rG/3)p(l + 3w), requiring only that 



1 Ref. has also addressed this question in the context of gravitiy with a time- varying speed 
of light. 
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w < —1/3. In (|1.2|) , as w approaches —1, the acceleration would disappear, com- 
ing into conflict with the observations. It must be emphasized however that since 
the connection between w and acceleration is no longer the same when eq. (|1 . 1|) is 
adopted, one should the modified expansion law. 

In the following we will present results for exponential quintessence potentials, 
since these proved to be the most promising for overcoming the horizon problem. 
As will be shown, the modified Friedmann equation enlarges the range of potential 
parameters which avoid a future horizon, but there is also an allowed range using the 
normal Friedmann equation. In the penultimate section it will be argued that these 
solutions do not require any more fine tuning than those arising from other models 
of quintessence. 



2. Quintessence evolution and event horizon 

The quintessence field Q with potential V has the usual equation of motion, Q + 
3HQ + ^ = 0. Its pressure and energy density are given by p Q = \Q 2 + V and 
Pq — \Q 2 — V. To explore cosmology both in the usual case and with the modified 
Friedmann equation ( |1 . 1|) we will introduce a parameter x = 0,1 and write the 
Hubble rate as 

H 2 = k\ ((1 + f )p r + p m + 1(1 + x)Q 2 + (1 - x)V) (2.1) 



where k 2 = 8irG/ (1+f) and p r and p m are the radiation and matter energy densities, 
respectively. It thus reduces to the standard equation when x = and the self-tuning 
one ( |1 . 1|) when x — 1. Guided by the CMB data, we assume the curvature term in 
H 2 to be absent. 

Rather than integrating with respect to time t, it is convenient to think in terms 
of redshift, 1 + z = l/a(t), where we take the present scale factor a(t ) to be unity. 
Further defining u = ln(l + z), and rescaling the fields and energy densities via 

Q = k x Q, pi = k x Hq pi, 

V = k 2 x Hu 2 V, H = H/H , (2.2) 

the quintessence equation of motion and Friedmann equation can be written in the 
dimensionless form 



Q" = H 



-2 



((1 + |)/5 r + \p m + (1 - X)V) Q'-(l- XQ' 2 ) -pr 



dQ\ 



(2.3) 



[l + l)Pr+p m + {l-x)V 



H 2 = y r " 1 \ ' , (2.4) 

l-i(l+x)Q' 2 

where primes denote 4-. The matter and radiation densities scale with u like p m = 
Pm,o£ 3u an d p r = p r ,oe 4u with respect to their present values at u = 0. These 
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dimensionless equations are well-suited to numerical integration, which is the main 
technique of our investigation. Notice that conventional time t has been eliminated, 
and u = ln(l + z) = — ln(a) now plays the role of the time variable. 

In comparing the properties of the quintessence field to observations, we will 
refer to the fractions of the critical energy density, and the equation of state. The 
former are defined as 



a = 


K xPi 


Pi . 


H 2 


H 2 ' 


% = 






Jp \ 


'}¥* 


which satisfy fl r 




+ n Q = 



(2.5) 

2 + (1 - X)V) = Q" + -1(1 - ( 2 -6) 

1. As for the quintessence equation of state, it is given 



by 



Q 2 -2V H 2 Q' 2 - 2V 

w = = — — 2.7 

Q 2 + 2V H 2 Q' 2 + 2V 



This follows from the fact that Q = —HQ'. 

In standard cosmology, w < —1/3 is the criterion for acceleration. But as noted 
in the introduction, the relation between w and acceleration is modified for the self- 
tuning scenario with x — 1. It is therefore useful to have another quantity indicative 
of acceleration, which can be more directly related to the observations of high-z 
supernovae. Let us first review what is actually constrained |]14||: it is the distance 



modulus (m — M = apparent minus absolute magnitude) of the SN versus its redshift, 
where m — M = 5 log 10 (c?L/Mpc) + 25, and the luminosity distance di(z) is given by 
(1 + z)H A(z), with 

. . , f*o dt f z dz' /•in(i+«) e « 



(z) a(t) Jo H(z') Jo H(u 

In a flat universe with only matter and cosmological constant components, Q m +Q/y = 
1, one would have 

Aa - r ,_ 3 „ d * (2.9) 



ii (x 3 n m + n A y/ 2 

The high-z SN results essentially try to measure A(z) as a function of z in order to 
fit We therefore define a phenomeno logical parameter, f2 A ff : 

Of = tt A such that A A (3/4) = A(3/4), (2.10) 

i.e., Q° A a is the amount of fl\ for pure cosmological constant that would be required 
to give the same value of A(z) as is produced by the quintessence model. This is a 
^-dependent definition, and we somewhat arbitrarily take the value z = 0.75 because 
this is roughly the redshift where the current SN data are which are the most sensitive 
to Q A . 
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The function A is also relevant for determining whether an event horizon appears 
in the spacetime, since the coordinate distance traveled by a photon between the 
present and a future (negative) value of z is precisely A(z). The criterion for an event 
horizon is that lira^-i A(z) < oo. In this case the photon travels a finite coordinate 
distance in an infinite time, and this determines the position of the horizon for the 
observer that emitted the photon: no signals originating from beyond that position 
will ever be able to reach him. For ease of representation, we will define the following 
measure of horizon formation: 

jj _ y din A f 0, a(t) ~ e Ht , horizon exists (2 11) 

a^oo ,9 In a | - — 1, a(t) ~ t q , no horizon if q < 1 

Here q = 3 ( 1 ^ w ) if the dominant component has equation of state p = wp. For 
example, D = 1/2 in a universe which behaves as though it is matter dominated 
(q = 2/3) at very late times. We will see below that D — 1/2 is the maximum value 
that arises in the quintessence models which we consider. 



3. Results 

We have examined some of the popular choices of the quintessence potential, V(Q), 
including inverse powers, exponentials, and combinations of the two. For avoiding the 



event horizon, the exponential potential (first considered in [15|, and subsequently 
in 



16[-||18|) seems most promising: 

V(Q) = V oe ~PQ (3.1) 

This potential has only a single free parameter, /3, 2 once the constraint that Q Q = 
1 — Q m is imposed for the present epoch, for this determines Vq. The statement is 
strictly true if one assumes that initially Q' — 0, since the initial value of Q itself 
can be absorbed into the definition of Vq. However, even if Q' ^ initially, the 



quintessence field converges to an attractor solution |]T6|-|jT8|, as shown in figure 
[I]: the two solutions corresponding to different initial conditions converge to the 
same functional form after some time. In this example, we numerically integrated 
the equation of motion for the case Q m = 0.25 starting from initial conditions at 
u = ln(l + z) — 12 (z — 1.6 x 10 5 ), during the radiation dominated era with initial 
conditions Q = 3 and Q' = or 1. Although the early behavior of Q is clearly 
affected by the difference in initial Q', both solutions join their common trajectory 
well before quintessence starts to dominate in the present era. Figure |2] shows the 
corresponding quintessence equation of state, w, in the two cases. In this example it 

2 in the notation of ref. |!q] , (3 — \/6a, and in that of refs. jl?], ^8[, (3 = V3X, assuming the 
normal Friedmann equation. 

3 The initial condition on Q is not significant since it can be absorbed into Vq. 
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is clear that quintessence can contribute to the acceleration during the period when 
w = — 1, whereas w — > —0.25 in the future, which is larger than —1/3 and therefore 
cannot cause acceleration. 




matter 

— - radiation 

— Quintessence, Q'=1 

— Quintessence, Q'=0 



10 5 -5 -10 -15 -20 

u = ln(1+z) 

Figure 1: £l r , £l m and Q Q as a function of u = ln(l + z) for the potential V = V^e 
{i.e., (3 = 2.6). The initial conditions at u = 12 are Q' = and Q' = 1, respectively. 
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10 5 -5 -10 

u = ln(1+z) 

Figure 2: The quintessence equation of state, w, as a function of u, for the same two sets 
of initial conditions as in figure [l]. 

Now let us consider the fate of the universe at very late times, for the interesting 
range of potential parameters. We find that a range of /3 values exists such that the 
cosmological expansion is accelerating today even though at late times it will revert 
to a power law, a ~ t q , with q < 2/3. This is illustrated by the solid curves of 
figure ^ which again were made assuming that fl m = 0.25 and the initial conditions 
Q — Q' — at u — 12. From the middle curve, the horizon formation parameter D 
(eq. ( [2. 11| )), we see that an event horizon is avoided if (3 £ 2.4. It is interesting to 



6 



note that for larger values of ft, D — > 1/2, just as it would for a matter-dominated 
universe, even though quintessence is dominating at late times. This happens because 
the solutions have the property that p Q ~ a~ 3 , just as though the universe was 
matter-dominated [IT] . 



However, very large values of ft do not yield sufficient acceleration at present 
times to be consistent with observations. This can be seen from the top and bottom 
curves, showing Q c f (defined in eq. ( |2.10|) ) and w (eq. (|2T7|)) respectively. The 99% 



confidence level SN limits are w < —0.5 for fl m = 0.25 and Q'f > 0.5 for a flat 
universe. It is difficult to apply the limit on w in a quantitative way in the present 
model because of the fact that w is changing very rapidly between the redshifts of 
z = 1/2, where much of the SN data is clustered, and the present, z = 0. This rapid 
variation can be seen in the example of figure § as well as by comparing the bottom 
sets of curves in figure |3|. On the other hand, the limit on f2^ ff can be applied in a 
straightforward way to give ft 2. 8. 4 Thus it is possible to satisfy the observational 
constraints without getting an event horizon if 2.4 <, ft < 2.8. 

We have also done the same analysis using the self-tuning Friedmann equation, 
shown by the dashed curves of figure || One sees that the range of ft for which no 
horizon forms is enlarged to ft £ 1-8. This agrees with the intuitive expectation that 
acceleration is reduced in this case, compared to the standard Friedmann equation. 
But at the same time, f2^ ff and \w\ are decreased, making it more difficult to obtain 
the observed acceleration. Since the connection between — and w is no longer the 
same as assumed in the analysis of the SN data, we again take advantage of the 
parameter. Demanding that > 0.5 gives ft < 2.4. The allowed range consistent 
with no horizon, 1.8 £ ft 2.4, is thus shifted and slightly widened relative to the 
normal case. 

It is easy to elucidate the origin of our loophole to the horizon-formation argu- 
ments put forward in references These analyses assume that quintessence is 
dominating very strongly in the present, or equivalently that its equation of state 
does not change from its present value. However this approximation is not valid for 
the solutions we have presented, such as in figure |l[ In the present epoch, u = 0, the 
equation of state w for these solutions is always going through a transition from — 1 
to some value greater than —1/3, so as to avoid the horizon. In the case of solutions 
with small ft < 2.4, it is also true that w changes near u = 0, but its final value at 
large times is w < — 1/3, leading to a horizon. 



4 Naively one might expect that 0^ = 1 — Vl m should be satisfied, which would rule out all values 
of j3 shown, but we remind the reader that fi^ ff is defined as the amount of real that would give 
the same amount of current acceleration as the given quintessence model. Since quintessence is less 
efficient at causing acceleration than is vacuum energy, it is not surprising that r2^ ff is less than fi Q . 
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Figure 3: From bottom to top: the quintessence equation of state w, at redshifts z = 1/2 
and at z = 0, the horizon formation parameter .D (eq. ( |2.11 )), and the SN effective Qa, 
fi^ (eq. (|2.10| )), as a function of the potential parameter /3. Solid lines are for the normal 
Friedmann equation {x = 0), dashed are for the self-tuning one (x = 1). 



4. Naturalness 



In this section we discuss the question of how much fine tuning is needed to obtain the 
desired solutions. The exponential potential has been somewhat disparaged because 
of the emphasis on solutions which reach the scaling regime, where w is constant, 
very early in the evolution. Such solutions are uninteresting in light of the current 
data because they maintain a constant value of Q Q . Since fl Q must be less than about 
15% at nucleosynthesis or during large scale structure formation, this would render 
its contribution too small to account for the present acceleration. But as pointed out 
in [19 1, this negative conclusion can be circumvented by assuming that quintessence 



is far from the late-time attractor solution in the not-too-distant past, so that w can 
evolve, which is exactly the situation for the solutions presented here. 

Does the fact that these quintessence solutions start out far from the late-time at- 
tractors make them less natural? We argue that this is not the case. All quintessence 
models require one tuning in order to achieve fl Q = 1 — Q m today, and this is the 
only one which we have invoked. The tuning is imposed as a particular value of the 



S 



combination 

V = V e-^, (4.1) 

where Qi is the initial value of Q. The value of Vo required to make Q Q = 1 — Q m 
today depends on the initial velocity, or equivalently Q\. In the very early universe 
the kinetic energy of quintessence typically dominates over its potential energy unless 
Qi is exactly zero (see eq. ( |2.6| )), so this amounts to a choice for the initial value of 
Q Qti = (1 + x)Qf/2. Figure |l] shows that whether f2 Qji = or 1, the recent evolution 
of the quintessence is identical, so long as V takes the right value. In this sense, 
we can say that the models under consideration are very insensitive to the initial 
conditions. This conclusion in no way depends on our choice for the initial time. 
Figure |] shows the evolution of for the same parameters as in fig. [l|, except now 
the initial redshift is taken to be 10 17 , corresponding to an initial temperature of 100 
TeV. Not only is the late-time evolution unaffected, but the choice of Vo is identical 
if Q\ = 0, and Vo only changes by a factor of 5 relative to the later initial condition 
if Q\ ± 0. 




40 20 -20 -40 



u = ln(1+z) = ln(1/a) 

Figure 4: Same as fig. [TJ, except starting at redshift z = 10 17 instead of z = 10 5 (and 
the additional initial condition Q[ = 0.04 is also shown, which corresponds to inflationary 
initial conditions with equipartition of energy). This demonstrates the insensitivity of the 
solution to the choice of the initial time. 

Another aspect of naturalness is the value of V required. If Q\ -C 1 (hence 
VL Q j <C 1), as would be natural in inflation if equipartition was realized |17|, [19[], then 
Vq must be of order the present critical density, Vo ~ (10~ 3 eV) 4 . (The evolution with 
such an initial condition, with fl Qti = 10~ 3 , hence Q\ = 0.04, is illustrated in figure 
^.) This looks unnaturally small in particle physics units, but one advantage of the 
exponential potential is that the smallness of Vo can be explained by a moderately 
large value of Qi. For example, if Vo ~ (TeV) 4 , then Qi should be of order 50M p , 
which is not such a disturbing hierarchy. 
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The final aspect of naturalness is how sensitive the solution is to small changes in 
Vq. Changing Vq causes a shift in the time when the quintessence and matter energy 
densities become comparable, which is related to the so-called coincidence problem: 
why is it that quintessence is just starting to dominate in the present epoch? It 
can be shown that scaling the potential by a factor of V new /Vo leads to the following 
dependence in the redshift of matter-quintessence equality, z m - q : 



The power 1/3 is just coming from the fact that Vq must be of order the critical 
density at z m - q , and density scales with redshift like (1 + z) 3 . Therefore, in some 
sense the value of z m _ q is rather insensitive to the value of Vo; the coincidence problem 
would be much worse if the dependence was through a higher power. 

5. Conclusion 

We have shown that it is possible to evade the cosmological event horizon which 
might pose a difficulty for deriving quintessence from string theory: for a range of 
2.4 <, p < 2.8 with the conventional Friedmann equation, or 1.8 <, (3 < 2.4 for the 
self-tuning variant, the exponential potential V(Q) = Voe~P KxQ gives this outcome. 
We also tried other kinds of potentials, such as inverse powers, V ~ Q~ p , but for 
these the development of a future horizon was found to be inevitable. 

We began this work with the idea that a self-tuning Friedmann equation might 
make it easier to avoid a cosmological horizon in a quintessential universe. The out- 
come is that self-tuning does not really make a big difference: eternal acceleration can 
be avoided with or without self-tuning. Of course self-tuning is still very interesting, 
because it allows us to work with a larger class of potentials, 



since only in the self-tuning case is the evolution completely insensitive to the value 
of V\. It is likely that there are other problems with self-tuning, since the strength 
of gravity on subgalactic scales is known to be consistent with the normal Planck 
mass, whereas gravity looks effectively weaker on cosmological scales in the self- 
tuning case. This is suggestive of the presence of an extra scalar component like a 
massive Brans-Dicke field, whose limited range accounts for the difference between 
the effective Planck mass at large and small distance scales. However the couplings 
of such a field to matter are very highly constrained by precision tests of general 
relativity in the solar system, like the precession of the perihelion of Mercury. 5 

5 1 thank Maxim Pospelov for discussions on this point. 




(4.2) 



V = V X + V e-P Q 



(5.1) 
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Regardless of self-tuning however, the class of solutions we have discussed seem 
sufficiently natural to warrant consideration as a strong candidate for the dark energy 
which is presently observed. It is expected that progress in the observations of high- z 
supernovae (the SNAP experiment |p20[) will soon be able to distinguish this kind 
of model from others through an accurate determination of the time dependence 
of the equation of state ||21||. The very large present time- dependence of w in the 



exponential models makes them particularly interesting in this respect. Figure |5] 
shows dw/dz at redshift z = 0.5 for the relevant range of (3. Thus if the SNAP 
experiment was to measure that dw/dz < —0.26 at z = 0.5, it would indicate that 
the universe will stop accelerating in the future and thus avoid an event horizon, in 
the context of the model discussed here. 



N 
"O 



-0.1 



-0.2 



-0.3 - 



-0.4 



normal Friedmann eq. 
self-tuning eq. 



-0.5 



1.6 1.7 1.8 1.9 



2.1 



2.2 2.3 2.4 2.5 2.6 



Figure 5: Rate of change of the quintessence equation of state, dw/dz, at z = 1/2 as a 
function of the potential parameter (3. In the region of (3 = 2.4, dw/dz ~ 0.23 — 0.2/3 for the 
normal Friedmann equation, and near = 1.8, dw/dz ~ —0.07 — 0.16/3 for the self-tuning 



one. 



Note added: as this work was being finished, ref. [22] appeared, which presents a 
different quintessence model that also avoids the future horizon. 
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